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ABSTRACT 

Optically-similar early-type galaxies are observed to have a large and poorly understood range in 
the amount of hot, X-ray-emitting gas they contain.To investigate the origin of this diversity, we 
studied the hot gas properties of all 42 early-type galaxies in the multiwavelength ATLAS^^ survey 
that have sufficiently deep Chandra X-ray observations. We related their hot gas properties to a 
number of internal and external physical quantities. To characterize the amount of hot gas relative 
to the stellar light, we use the ratio of the gaseous X-ray luminosity to the stellar AT-band luminosity, 
we also use the deviations of from the best-fit Lx^^^-Lk relation (denoted ALxg„)- 

We quantitatively confirm previous suggestions that various effects conspire to produce the large 
scatter in the observed Lx!Lx relation. In particular, we find that the deviations ALx^^^ are most 
strongly positively correlated with the (low rates of) star formation and the hot gas temperatures in 
the sample galaxies. This suggests that mild stellar feedback may energize the gas without pushing 
it out of the host galaxies. We also find that galaxies in high galaxy density environments tend to 
be massive slow-rotators, while galaxies in low galaxy density environments tend to be low mass, 
fast-rotators. Moreover, cold gas in clusters and fields may have different origins. The star formation 
rate increases with cold gas mass for field galaxies but it appears to be uncorrelated with cold gas for 
cluster galaxies. 

Subject headings: X-rays: galaxies: luminosity - galaxies: ISM - galaxies: elliptical and lenticular 
Clusters of galaxies: intracluster medium 


1. INTRODUCTION 

Present day early-type galaxies (elliptical and lentic¬ 
ular galaxies; ETCs) tend to contain substantial atmo¬ 
spheres of hot gas. The thermal and chemical properties 
of this gas reservoir preserve much of the history of the 
galaxy. Such hot gas can also regulate star formation and 
impact the evolution of the host galaxy. In particular, 
there have been a growing number of studies that show 
ETCs may contain sizable cold gas halos (Oosterloo et 
al. 2010; Young et al. 2011; Serra et al. 2012). Lagos 
et al. (2014) found through numerical simulations that 
90% of the cold gas in ETCs may be supplied by the ra¬ 
diative cooling of their hot atmospheres. The OVH line 
has been detected in the ISM of several ETCs, revealing 
the existence of weak cooling flows (Pinto et al. 2014). 
Given that hot gas is heavily involved in the formation 
and evolution of ETCs, it is crucial to understand what 
mechanisms are capable of modifying the hot gas con¬ 
tent. 

The primary source of the hot ISM is gas lost from 
aging stars, while at least some fraction comes from mi¬ 
nor mergers and accretion of the intergalactic medium 
(Mathews & Brighenti 2003). Such gas is heated through 
supernova explosions and the thermalization of stel¬ 
lar motions. It radiates mainly in X-rays via thermal 
bremsstrahlung and line emission. Galaxies with compa¬ 
rable stellar masses should have similar amounts of gas 
deposited into the ISM. Thus, we expect to observe a 
tight relation between the X-ray luminosity (Lx) and 
the optical luminosity (Lopt)- In contrast, a puzzlingly 
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large variation in X-ray luminosity among ETCs of simi¬ 
lar optical luminosities has been reported in observations. 
The discussion of this discrepancy dates back to the time 
of the Einstein Observatory (e.g. Canizares et al. 1987; 
Fabbiano et al. 1992). As one of the early investigations 
using the modest ROSAT measurements, O’Sullivan et 
al. (2001) found a scatter in Lx/Aopt of up to two orders 
of magnitude for 430 nearby ETCs. Later, the Chan¬ 
dra X-ray Observatory, with its high sensitivity, large 
spectral range, and in particular its superb 0.5" angular 
resolution, found that ETGs contain detectable X-ray 
emitting point sources such as low mass X-ray binaries 
(LMXBs), in addition to hot gas (e.g. Sarazin 2000). Re¬ 
cent studies show that additional stellar components such 
as cataclysmic variables and active binaries (CV/ABs) 
contribute to the diffuse X-ray emission as well (e.g. 
Revnivtsev et al. 2008). Galaxies with low Lx/Aopt 
generally have a larger fraction of their flux provided 
by stellar point sources than do X-ray luminous gas-rich 
galaxies. After excluding the contamination of LMXBs 
and other stellar emission, current studies with Chandra 
observations reveal an even larger scatter in the gaseous 
content relative to stars, as measured by the Axg^^^/Aopt 
relation for ETGs (e.g. Boroson et al. 2011). We might 
naively expect that the Ax^^^/Aopt relation may be re¬ 
lated to galaxy ages, since older galaxies have more time 
to accumulate gaseous atmospheres; however the scatter 
in Ax /Aopt greatly exceeds the range in ages for these 
ETGs.“ 

Energy feedback may modify the hot gas content of 
ETGs, since gas can be driven from galaxies by super¬ 
novae (SNe) heating and active galactic nuclear (AGN) 
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activity. David et al. (2006) found that heating by Type 
la supernova (SNIa) is energetically sufficient to gener¬ 
ate galactic winds in these galaxies even if the present 
SNIa rate is overestimated. AGN heating may also in¬ 
stigate gaseous outflows. X-ray cavities filled with radio 
emission have been frequently detected in the ISM of el¬ 
liptical galaxies (e.g. Finoguenov & Jones 2002, Forman 
et al. 2005), which indicates the effects of AGN activity. 
If the feedback is relatively mild, gas may just be re¬ 
distributed out to larger radii (and eventually fall back) 
and the feedback energy could heat up more gas into an 
X-ray emitting phase (Chevalier & Clegg 1985). These 
internal feedback processes are particularly important for 
isolated galaxies, since environmental effects such as ram 
pressure stripping and pressure confinement by intraclus¬ 
ter gas would be insignificant. 

Deeply connected to the feedback scenario is the de¬ 
pendence of a galaxy’s gaseous content upon the depth 
of the galaxy’s gravitational potential: it is easier for 
more massive galaxies to retain their gas against loss 
mechanisms such as ram pressure stripping and inter¬ 
nally driven outflows. Mathews et al. (2006) studied a 
number of massive ETGs and found that their hydro¬ 
static masses seem to correlate with Lx /Topt; they con¬ 
cluded that the role of feedback becomes more impor¬ 
tant in less massive galaxies. However, the total masses 
of ETGs can be very difficult to measure accurately. In 
some cases, hot gas temperatures and stellar velocity dis¬ 
persions can be regarded as proxies for total masses. The 
depths of gravitational potentials can also be influenced 
by the rotation and flattening of galaxies. It has been no¬ 
ticed since the Einstein Observatory epoch that flattened 
ETGs have smaller Lx /Topt compared to their rounder 
counterparts (Eskridge et al. 1995). Sarzi et al. (2013) 
demonstrated that slowly rotating galaxies have larger 
X-ray halos, although it has always been difficult to dis¬ 
entangle rotation and flatness, since fast rotational sys¬ 
tems would also be flatter due to centrifugal forces. The 
ATLAS^'^ project proposes a kinematic classification of 
ETGs into fast- and slow-rotators to distinguish “true 
ellipticals” from “true lenticulars;” this kinematic clas¬ 
sification is less subject to the projection effects which 
confounded prior morphological studies (Emsellem et al. 
2011 ). 

Environmental influences are another widely discussed 
possible cause of the scatter in Lxi^^^/Lopt (White & 
Sarazin 1991; Mathews & Brighenti 1998; Brown & Breg- 
man 2000). Most galaxies, in particular early-type, re¬ 
side in groups and clusters. The hot gas content of such 
galaxies can be affected by their environments through 
various processes, including interactions between galax¬ 
ies (harassment), interactions between galaxies and intr¬ 
acluster gas (ram pressure stripping and pressure confine¬ 
ment), and interactions between galaxies and the gravita¬ 
tional potential of the cluster (tidal stripping). Perhaps 
the most commonly discussed case is ram pressure strip¬ 
ping, the process through which the hydrodynamic drag 
of intracluster gas pulls the ISM out of galaxies mov¬ 
ing through it (Gunn & Gott 1972). Simulations suggest 
that nearly all cluster galaxies experience ram pressure 
stripping throughout their lifetime (Bruggen & De Lucia 
2008). Indeed, nothing is more convincing than witness¬ 
ing a galaxy losing its hot gas. There are several X-ray 
observations of on-going ram pressure stripping in ETGs. 



Fig. 1.— K-band luminosity distribution of all 260 galaxies in 
ATLAS^^ and 42 galaxies in our sample. K-band luminosities 
shown in this figure are taken from Cappellari et al. (2011). 

These observations are characterized by hot gas tails dis¬ 
placed from, and trailing behind, the host galaxy’s stellar 
distributions (M86 - Randall et al. 2008; NGC 4472 - Ir¬ 
win & Sarazin 1996; NGG 1404 - Machacek et al. 2005; 
NGG 1400 - Su et al. 2014). Most such studies focus only 
on individual galaxies. The next logical step is to under¬ 
take a statistical study of a sample of ETGs experiencing 
varying amounts of ram pressure. 

Our primary goal is to understand what causes the 
large scatter in the /-^opt relation, using a statisti¬ 
cally large sample of ETGs. In this paper we investigate 
the hot gas content of all 42 ETGs in the ATLAS^° sur¬ 
vey that have sufficiently deep Chandra X-ray observa¬ 
tions. We look for correlations between Lx^^^/Lk and 
various internal and external factors which may affect 
the hot gas content of ETGs, in order to infer their rel¬ 
ative importance. In addition, 15 of the 42 galaxies in 
this sample reside in the Virgo Gluster, which allows us 
to statistically study the effects of ram pressure stripping 
in a cluster environment. 

We assume Hq = 70 km s“^ Mpc“^, = 0.7, and 
Dm = 0.3. Throughout this paper uncertainties are 
given at the Icr confidence level unless otherwise stated. 
We report our sample selection process and introduce 
their properties in §2; observations and data reduction 
are described in §3; we report our results in §4; we dis¬ 
cuss the implications of our results in §5; we summarize 
our main conclusions in §6. In the Appendix, we exam¬ 
ine statistical and systematic uncertainties and provide 
additional figures. 


2. SAMPLE SELECTIONS 

We want to study a statistically significant sample 
of ETGs that is as diverse as possible and with pub¬ 
lished values of multiple relevant physical characteristics 
of the galaxies. Our main goal is to relate the hot gas 
properties of these galaxies to those observables which 
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can reflect their formation and evolution from differ¬ 
ent aspects. ATLAS^^ is a volume-limited (1.16x10^ 
Mpc^), multi-wavelength survey of a complete sample 
{Lk > 6 X 10® Lq) of ETGs (Cappellari et al. 2013). 
To our knowledge, it provides the richest and most com¬ 
plete information on the global stellar kinematics, dy¬ 
namics, and different phases of the interstellar medium 
(but not X-ray emitting hot gas) of a large sample of 
ETGs. Even more importantly, the values of these ob¬ 
servables were derived through a uniform analysis. We 
aim to study all galaxies in ATLAS^° that have Chandra 
observations that are deep enough to constrain their hot 
gas properties. 

We require that each galaxy has been observed 
with Chandra for at least 15 ksecQ Our screen¬ 
ing criteria yielded 42 ETGs, as listed in Table 1. 
We take their distances and half-light radii {r^) from 
ATLAS^'® (Cappellari et al. 2013). All galaxies in our 
sample reside at distances between 10 and 31 Mpc, with 
their ranging from 0.5 to 8.6 kpc. Galaxies in our sam¬ 
ple are diverse in luminosity and environment and are 
representative of ETGs. From the literature we obtain 
stellar ages, stellar velocity dispersions (tr), atomic gas 
masses (Mhi), molecular gas masses (Mna), total masses 
(within Ire) (Aftot), a rotational parameter (A), and 
galaxy ellipticities (e), with most values from ATLAS®^ 
and measured with consistent analyses. We take the star 
formation rate (SFR) of almost all galaxies in our sample 
from Amblard et al. (2014). These values were derived 
from their spectral energy distributions (SEDs), fitting 
from the UV to millimeter wavelengths. If not listed in 
Amblard et al. (2014), we take SFR values from Davis 
et al. (2014), which are derived from WISE 22/im data. 
We list these observables in Tables 2 and 3. 

It should be noted that Chandra archived galaxies do 
not form a complete sample. In Figure [l] we compare 
the K-band luminosity distribution of 42 galaxies in our 
sample and 260 galaxies in the ATLAS®^ survey. Very 
faint galaxies {Lk < 10^® Li^©) have been underrepre¬ 
sented in our sample. To our knowledge, this is however 
the largest sample available that is suitable for associat¬ 
ing hot gas properties of ETGs with internal and external 
factors that may be related to galaxy evolution. 

In order to probe how different factors affect the 
ISM of different galaxies, we also subdivide our sam¬ 
ple in three ways, based on the E/SO dichotomy {X/y/e), 
galaxy mass, and environment. In the E/SO subdivi¬ 
sion, we use the stellar dynamical criterion found by the 
ATLAS®'® project to distinguish “true elliptical” galax¬ 
ies from “true lenticular” galaxies in our sample. The 
ATLAS®'® project found that a combination of a rota¬ 
tional parameter (A) and the galaxy ellipticity (e) could 
be used to distinguish slow rotators (“true ellipticals,” 
which have Xjy/e < 0.31) from fast rotators (“true lentic- 
ulars,” which have X/y/e > 0.31). Using this criterion for 
our sample, we will distinguish 12 “true elliptical” galax¬ 
ies from 30 “true lenticular” galaxies. We also divide our 
sample into high mass and low mass systems: we re- 

^ We did not include NGC 4486A since this galaxy is too close 
to M87. It is extremely difficult to disentangle its X-ray emission 
from that of M87. We did not include M87 because it is the central 
galaxy of the Virgo cluster, so it is difficult to separate its galactic 
gas from intracluster gas. 


gard the 20 galaxies in our sample with Mtot > 10" Mq 
as high mass galaxies, while the remaining 22 galaxies 
are regarded as low mass. Finally, we divide our sample 
environmentally: galaxies with 15 or fewer SDSS galax¬ 
ies nearby (Column 4 in Table 3; see §3.5) are classi¬ 
fied as field galaxies, while galaxies with more neighbors 
are classified as belonging to groups and clusters. Four 
galaxies in our sample do not have good SDSS coverage, 
so we determine their environments differently: accord¬ 
ing to the literature, we consider NGC 1023 as a galaxy 
in a group and we regard NGC 821, NGC 1266, and 
NGC 7457 as being in relatively isolated environments 
(Chernin et al. 2010; Arnold et al. 2014; Malchaney & 
Jeltema 2010; Alatalo et al. 2014). Consequently, 21 
galaxies in our sample are classified to be in high galaxy 
density environments and the other 21 galaxies are clas¬ 
sified to be in low galaxy density environments. 


3. OBSERVATIONS AND DATA REDUCTION 

Chandra is the best instrument to study the hot gas 
properties of ETGs, which have kT < 1 keV. Its superb 
spatial resolution can resolve out LMXBs and central 
active galactic nuclei. The emission of other stellar com¬ 
ponents are also best calibrated with Chandra. We used 
CIAO 4.5 and CALEB 4.5.8 to reduce ACIS-S and ACIS- 
I data. All data were reprocessed from level 1 events, so 
that the latest, consistent calibrations were used. Only 
events with grades 0, 2, 3, 4, and 6 are included. We also 
removed bad pixels, bad columns, and node boundaries. 
We filtered background flares with the light curve filter¬ 
ing script lc_clean. The effective exposure times are 
shown in Table 1. Point sources were detected in a 0.3- 
7.0 keV image with wavdetect, supplied with a 0.5 keV 
exposure map. The detection threshold was set to 10“® 
and the scales of wavdetect ranged from 1 to 8, in steps 
increasing by a factor of Point sources, including 
galactic nuclei, were removed. 

3.1. Regions and Background 

We adopt the optical effective radii re for the sample 
galaxies as listed by ATLAS®'®, which were taken from 
the Third Reference Catalogue of Bright Galaxies (RC3, 
de Vaucouleurs et al. 1991). We chose an X-ray extrac¬ 
tion aperture of 2 re to determine the X-ray luminosity 
(Lxgas) and mass {Mx^^J of the hot gas. Local back¬ 
ground, extracted from a region away from the source 
region on the same CCD chip, was used in the X-ray 
spectral analysis of most galaxies. The area of the lo¬ 
cal background was chosen to be at least twice the area 
of the source region to ensure a sufficient S/N ratio for 
background subtraction. In some cases the ISM emission 
fills the entire S3 chip. We adopt other methods to esti¬ 
mate the background emission for galaxies with effective 
radii larger than 0.9'. For galaxies also observed with 
the SI chip, we used “stowed background” dat£[^for the 
spectral fitting. We fit a spectrum extracted from the SI 
chip with a stowed background of the same region on the 
SI chip to determine the surface brightness of cosmic X- 
ray and Galactic emission background, since the SI chip 

® particle background inside the detector observed with ACIS 
stowed, http://cxc.harvard.edu/contrib/maxim/stowed 
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TABLE 1 

Chandra Observational Log 


Name 

Instrument 

Obs-ID 

Exposure (ks) 

IC1024 

ACIS-S 

14901 

19.8 

NGC 821 

ACIS-S 

4006, 4408, 5691,5692, 6310, 6313, 6314 

13.9,13.6,39.6, 27.6, 31.8,49, 39.6 

NGC 1023 

ACIS-S 

8197,8198,8464,8465 

46.8,47.3,42.4,44.2 

NGC 1266 

ACIS-S 

11578 

29.5 

NGC 2768 

ACIS-S 

9528 

61.8 

NGC 2778 

ACIS-S 

11777 

29.2 

NGC 3377 

ACIS-S 

2934 

38.8 

NGC 3379 

ACIS-S 

1587, 7073, 7074, 7075, 7076 

29.9, 81.6, 67.9, 79.4, 66.5 

NGC 3384 

ACIS-S 

13829 

92.5 

NGC 3599 

ACIS-S 

9566 

19.7 

NGC 3607 

ACIS-I 

2038 

27 

NGC 3608 

ACIS-I 

2038 

27 

NGC 3665 

ACIS-S 

3222 

17.9 

NGC 4203 

ACIS-S 

10535 

41.6 

NGC 4261 

ACIS-S 

9569 

99 

NGC 4278 

ACIS-S 

4741, 7077, 7079, 7080, 7081, 7808 

35,111, 106,56,112, 52 

NGC 4342 

ACIS-S 

4687, 12955 

35.4, 68.1 

NGC 4365 

ACIS-S 

2015, 5921, 5922, 5923, 5924, 7224 

39, 38, 39, 35, 25, 10 

NGC 4374 

ACIS-S 

5908, 6131, 803 

45.3, 38.9, 27.3 

NGC 4382 

ACIS-S 

2016 

39.7 

NGC 4406 

ACIS-S 

318, 963 

11, 14 

NGC 4458 

ACIS-S 

14905 

29.4 

NGC 4459 

ACIS-S 

2927, 11784 

9.8, 29.8 

NGC 4472 

ACIS-S 

321, 11247 

31.5, 38.6 

NGC 4473 

ACIS-S 

4688 

25.5 

NGC 4477 

ACIS-S 

9527 

35.0 

NGC 4494 

ACIS-S 

2079 

21.7 

NGC 4526 

ACIS-S 

3925 

39.6 

NGC 4552 

ACIS-S 

2072 

51.3 

NGC 4564 

ACIS-S 

4008 

15.7 

NGC 4596 

ACIS-S 

11785 

31.0 

NGC 4621 

ACIS-S 

2068 

21.3 

NGC 4636 

ACIS-S 

323 

43.5 

NGC 4649 

ACIS-S 

785, 8182, 8507 

15.8, 36.9, 17.1 

NGC 4697 

ACIS-S 

784, 4727, 4728, 4729, 4730 

38, 40, 34, 22, 35 

NGC 4710 

ACIS-S 

9512 

28.8 

NGC 5422 

ACIS-S 

9511, 9772 

17.9, 18.4 

NGC 5576 

ACIS-S 

11781 

29.7 

NGC 5813 

ACIS-S 

12952 

143 

NGC 5846 

ACIS-S 

788, 7062, 8448, 8449 

21,23, 8, 19 


ACIS-I 

7923 

87 

NGC 5866 

ACIS-S 

2879 

30.6 

NGC 7457 

ACIS-S 

11786 

28.7 


is more offset and less contaminated by source emission. 
We then fit the spectrum extracted from the S3 chip with 
a corresponding stowed background by adding scaled X- 
ray background components obtained with the SI chip 
to the fitting. If the SI observation was not available, we 
used ACIS “blank-sky” datgj^for a background estimate. 
Since the amplitude of the instrumental background com¬ 
ponents vary with time, we scale the background count 
rate using the 10-12 keV count rate ratio of the source 
and background data. 

3.2. Spectral analysis 

We grouped spectra to have at least one count per 
energy bin and adopted C-statistics for all our spectral 
analyses. Redistribution matrix files (RMFs) and ancil¬ 
lary response files (ARFs) were generated for each re- 

http://cxc.harvard.edu/contrlb/maxim/aclsbg/ 


gion using the specextract tool. Energy bands were 
restricted to 0.7-7.5 keV, where the responses are best 
calibrated. We performed spectral analysis with Xspec 
12.7.0. The model we adopted to fit the diffuse emission 
in each galaxy is phabs * (apec -f powerlaw -L mekal -|- 
powerlaw). The absorbing column density Nh, asso¬ 
ciated with the photoelectric absorption model compo¬ 
nent phabs, was fixed at the Galactic value in the line 
of sight to each galaxy (Dickey & Lockman 1990). The 
apec component represents thermal emission from the 
hot gas. The first powerlaw component, with an index 
fixed at 1.6, represents the contribution from unresolved 
LMXBs (Irwin et al. 2003). In addition to hot gas and 
unresolved LMXBs, faint stellar X-ray sources such as 
cataclysmic variables (CVs) and coronally active binaries 
(ABs) also contribute to the X-ray flux. Revnivtsev et 
al. (2007, 2008, 2009) calibrated the X-ray emission from 
such old stellar populations in several extremely gas-poor 
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TABLE 2 

Properties of early-type galaxies in our sample from ATLAS3D 


Name 

Distance^ 

(Mpc) 

' e 

(kpc, arcmin) 

Mtot*’ 

(lOiOM©) 

Rotation 

Ellipticity 

A/v^ 

Mhi'* 

(10«M©) 

(1O«M0) 

(7^ 

(km s~^) 

IC1024 

24.2 

1.293, 0.187 

1.469 

0.691 

0.587 

0.902 

- 

3.981 

77.98 

NGC 821 

23.4 

4.507, 0.664 

12.455 

0.273 

0.392 

0.436 

<0.081 

<0.331 

179.5 

NGC 1023 

11.1 

2.575, 0.798 

6.592 

0.391 

0.363 

0.649 

19.498 

0.062 

166.7 

NGC 1266 

29.9 

2.891, 0.340 

2.576 

0.638 

0.193 

1.452 

- 

19.498 

79.07 

NGC 2768 

21.8 

6.499, 1.052 

34.198 

0.253 

0.472 

0.368 

0.646 

0.437 

198.2 

NGC 2778 

22.3 

1.696, 0.264 

3.141 

0.572 

0.224 

1.209 

<0.115 

<0.302 

132.1 

NGC 3377 

10.9 

1.870, 0.591 

2.938 

0.522 

0.503 

0.736 

<0.033 

<0.091 

128.2 

NGC 3379 

10.3 

1.975, 0.664 

8.222 

0.157 

0.104 

0.487 

<0.031 

<0.052 

185.8 

NGC 3384 

11.3 

1.773, 0.539 

3.664 

0.407 

0.065 

1.596 

0.178 

0.129 

130.0 

NGC 3599 

19.8 

2.236, 0.391 

0.989 

0.282 

0.08 

0.997 

<0.107 

0.234 

63.68 

NGC 3607 

22.2 

4.163, 0.648 

21.979 

0.209 

0.185 

0.486 

<0.083 

2.630 

206.5 

NGC 3608 

22.3 

3.158, 0.492 

9.036 

0.043 

0.19 

0.099 

0.145 

<0.380 

169.0 

NGC 3665 

33.1 

4.876, 0.515 

35.975 

0.41 

0.216 

0.882 

<0.269 

8.128 

216.3 

NGC 4203 

14.7 

2.095, 0.492 

4.018 

0.305 

0.154 

0.777 

14.125 

0.245 

129.1 

NGC 4261 

30.8 

5.539, 0.634 

52.723 

0.085 

0.222 

0.18 

0.5^ 

<0.479 

265.5 

NGC 4278 

15.6 

2.350, 0.527 

11.912 

0.178 

0.103 

0.555 

6.310 

<0.282 

212.8 

NGC 4342 

16.5 

0.518, 0.110 

3.319 

0.528 

0.442 

0.794 

- 

<0.174 

242.1 

NGC 4365 

23.3 

5.831, 0.875 

33.497 

0.088 

0.254 

0.175 

<0.4f 

<0.417 

221.3 

NGC 4374 

18.5 

4.653, 0.875 

38.459 

0.024 

0.147 

0.063 

<0.182 

<0.170 

258.2 

NGC 4382 

17.9 

5.656, 1.101 

28.054 

0.163 

0.202 

0.363 

<0.093 

<0.229 

179.1 

NGC 4406 

16.8 

7.512, 1.555 

39.811 

0.052 

0.211 

0.113 

1 

<0.251 

190.5 

NGC 4458 

16.4 

1.838, 0.391 

1.076 

0.072 

0.121 

0.207 

<0.081 

<0.204 

88.51 

NGC 4459 

16.1 

2.817, 0.605 

8.299 

0.438 

0.148 

1.139 

<0.081 

1.738 

158.1 

NGC 4472 

17.1 

7.882, 1.592 

59.566 

0.077 

0.172 

0.186 

0.5^ 

<0.178 

250.0 

NGC 4473 

15.3 

2.000, 0.449 

8.472 

0.229 

0.421 

0.353 

<0.072 

<0.117 

186.6 

NGC 4477 

16.5 

3.050, 0.648 

8.770 

0.446 

0.135 

1.214 

<0.089 

0.347 

148.9 

NGC 4494 

16.6 

3.943, 0.816 

9.840 

0.212 

0.173 

0.51 

<0.069 

<0.178 

150.0 

NGC 4526 

16.4 

3.551, 0.744 

17.498 

0.453 

0.361 

0.754 

<0.198 

3.890 

208.9 

NGC 4552 

15.8 

2.517, 0.565 

15.922 

0.049 

0.047 

0.226 

<0.074 

0.191 

224.4 

NGC 4564 

15.8 

1.517, 0.340 

3.828 

0.619 

0.56 

0.827 

<0.081 

0.178 

154.5 

NGC 4596 

16.5 

3.050, 0.648 

8.204 

0.639 

0.254 

1.268 

<0.135 

0.204 

125.6 

NGC 4621 

14.9 

3.080, 0.711 

12.882 

0.291 

0.365 

0.482 

<0.072 

<0.135 

197.7 

NGC 4636 

14.3 

6.167, 1.485 

24.889 

0.036 

0.094 

0.117 

8.lf 

<0.074 

181.6 

NGC 4649 

17.3 

5.453, 1.101 

52.360 

0.127 

0.156 

0.322 

<0.151 

<0.275 

267.9 

NGC 4697 

11.4 

3.379, 1.028 

11.695 

0.322 

0.447 

0.482 

<0.6^ 

<0.072 

169.4 

NGC 4710 

16.5 

2.368, 0.503 

5.768 

0.652 

0.699 

0.78 

0.069 

5.248 

104.7 

NGC 5422 

30.8 

3.115, 0.356 

8.892 

0.6 

0.604 

0.772 

0.741 

<0.603 

157.4 

NGC 5576 

24.8 

2.609, 0.365 

7.568 

0.102 

0.306 

0.184 

<0.4f 

<0.398 

155.2 

NGC 5813 

31.3 

8.614, 0.959 

39.08 

0.071 

0.17 

0.172 

<i.of 

<0.490 

210.9 

NGC 5846 

24.2 

6.783, 0.981 

36.559 

0.032 

0.062 

0.129 

3.5 

<0.603 

223.4 

NGC 5866 

14.9 

2.621, 0.605 

10.023 

0.319 

0.566 

0.424 

0.091 

2.951 

157.0 

NGC 7457 

12.9 

2.251, 0.605 

1.652 

0.519 

0.47 

0.757 

<0.041 

0.091 

74.64 


Note. — (a). Cappellari et al. (2011); (b). Total mass measured within IPe (Cappellari et al. 2013); (c). Emsellem et al. (2011); (d). 
Serra et al. (2012); (e). Young et al. (2011); (f). Serra & Oosterloo (2010); (g). Lucero & Young (2013). 


galaxies. Revnivtsev et al. (2008) found a nearly univer¬ 
sal relation for the unresolved X-ray emissivity per Lk 
in old stellar populations: LxILk = 5.9 x lO^^ergs”^ 
Lkq~^, where Lx is in the 0.5-2.0 keV band. We 
used this relation to derive X-ray estimates for these 
stellar sources from the galaxies’ RT-band luminosities 
Lk within 2re- The mekal+powerlaw components rep¬ 
resent these CV/ABs sources, where the mekal ther¬ 
mal emission temperature is fixed at 0.5 keV (and its 
abundance is fixed at solar), and the power law index is 
fixed at 1.9 (Revnivtsev et al. 2008). The ratio of the 
fluxes of the mekal and powerlaw components was set 
to 2.03 (Revnivtsev et al. 2008). In our spectral analy¬ 
sis, we hxed the CV/ABs components at the estimated 
flux based on the Lk of each galaxy. The redshift of 
each galaxy was taken from NED and fixed. We let the 


gaseous apec model temperature, metallicity, and nor¬ 
malization free to vary. For those galaxies which did 
not have well-constrained metallicities, we fixed them at 
0.7 Zq, the average of the hot gas metallicities of ETCs 
{kT < IkeV) studied by Su & Irwin (2013). 

3.3. Determination of Txg„ and 

The X-ray luminosities estimated in this paper 

are hot gas luminosities in the 0.1 —2.0 keV energy range. 
X-ray emission from CV/ABs and unresolved LMXBs 
are explicitly excluded, having been removed spectrally 
as outlined above. We list the best-ht Lx^^^ and gas 
temperature Tx for each galaxy in Table 3. Assuming 
a spherical distribution for the hot gas, we obtained the 
volume of the extraction region for each galaxy. We de¬ 
rived the average hot gas density from the best ht nor- 
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TABLE 3 

Other properties of early-type galaxies in our sample and their spectral fitting results 


Name 

Morphology* 

Age 

SDSS 

Virgo^ 

SFR'> 

Lx * 

Lk^ 

Tx* 




(Gyr) 



(Me/yr) 

(10"° ergs *) 

(10*°L^b) 

(keV) 

(10°Mq) 

IC1024 

so 

- 

5 

0 

U.OZ^_o 154 

u.uyo_o 032 

1.116 

- 

0.112+™;;- 

NGC 821 

E6 

7.5^ 

- 

0 

< 0.013'=! 

u.uuoy_o 0039 

8.930 

0 1 k + U.85* 
U- 10_o 05 

0.198l°;“; 

NGC 1023 

SBO 

4.7® 

- 

0 

U.UiO_Q QQg 

u.ooy_o 348 

6.550 

0 -1 K+0.01 

U--LO_o 02 

1 003+0-450 
l.UU9_o 314 

NGC 1266 

SBO 

l.S*’ 

- 

0 

1 “^1 q+0.920 

i.OiO_Q 542 

n 988+0-071 
U.2«»_o 049 

3.002 

n Qy+O-OO 

y-o<^_0.09 

u.oyz_o 062 

NGC 2768 

E6 

11.2™ 

5 

0 

u.uoo_Q 014 

i.z^y_0 183 

15.600 

0 37+0-04 

4 907+0-462 
^.yu/_o 374 

NGC 2778 

E 

5.4“" 

4 

0 

U.UU8_o 004 

0 007+0-010 
U.UU/ _o 007 

2.022 

- 

0 045+0-026 
U.U40_o 045 

NGC 3377 

E5-6 

3.7“^ 

12 

0 

U.UU4_0_oo3 

U.UU9Z_o 0032 

2.274 

^•^^-0.07 

0 041+0-030 

U.U4i_0 041 

NGC 3379 

El 

8.6^ 

12 

0 


0 01 0+0-003 
U.UiU_o 0025 

5.777 

0 9K+0.03* 

U-^0_o 02 

U.UOD_o 009 

NGC 3384 

SBO 

3.2'^ 

12 

0 

0.008t°;°w 

0 007+0-004 

U.UU/ _o 004 

4.038 

0 9K+0.30 

0 071+0-018 
u.u/ i_o 021 

NGC 3599 

SAO 

- 

13 

0 

0 Q7Q+0 021p 

U.U/O_0 016 

0 01 c+O.OlO 
U.UiO_o 007 

1.575 

- 

0.095l°;“f« 

NGC 3607 

SAO 

3.6'^ 

13 

0 

u.u^z_o 018 

0 740+0-106 
u. ^^O_o 084 

12.380 

0 45+0-03 
U-4O_0 05 

2.078t°;;« 

NGC 3608 

E2 

6.9*^ 

13 

0 

0 m c;+o.oo7 
U.U10_o_oo5 

0 3C0+O.O75 

U.900_o_o60 

4.793 

Q oq+0.07 

0 995+0-100 

u.yyo_o 087 

NGC 3665 

SAO 

9.9’^ 

1 

0 

U.iuy_o 032 

1 919+0-385 

i.yiy_0 287 

18.183 

0.45l°;“ 

4 347+0-416 
-0.338 

NGC 4203 

SABO 

- 

2 

0 

0 Q47+0 035 

U.U4t/ _o 020 

0 039+0 014 
u.uoy_o 014 

3.914 

o.83i°-;° 

0 1 C0+O.O25 
U.iOD_o 030 

NGC 4261 

E2-3 

15.5^ 

51 

0 

Q . .-.+0.064 
^ -0.042 

4-26ll°;°“ 

21 

0 81 +0-01 
d-ol_o 01 

8 1 qj?+0.036 

5.i»b_0 036 

NGC 4278 

El-2 

10.7'^ 

17 

0 

u.u±y_o 007 

0.0881°;“^ 

6.3 

0 65+0-06 

U-O0_o 06 

0 329+0-013 

u.ozy_o 023 

NGC 4342 

SO 

- 

32 

1 

- 

0 0[5c:+O.OO7 

U.UOO_o_oo6 

1.25 

u-oz_o 01 

0 067+0-004 
U.UO/_Q QQ4 

NGC 4365 

E3 

5.9'i 

43 

0 

U.UOU_o 017 

0 544+0-041 

U.O^^_0_o40 

16 

0 59+0-02 
u-oy_o 02 

2 447+0-090 

-0.091 

NGC 4374 

El 

12.2^ 

46 

1 

0 f)[5fi+0 026 

U.UO«_o_oi8 

5 423+0-344 

o.^zo_o 531 

19.76 

0 76+0-004 

O'- '^-0.004 

6 044+0-296 
o.u^^_0 303 

NGC 4382 

SAO 

1.6'= 

15 

1 

U.UUZ_o 002 

1 o'7q+ 0.251 
i.0/0_o 259 

20.82 

0 OO+0.02 

y-'JO_0.02 

c 0^1+0.467 

0.90i_Q 53Q 

NGC 4406 

E3 

11® 

28 

1 

- 

9.9881;-“; 

26.82 

0 81 +0-01 
d-ol_o.oi 

14 81 +0.921 

14.5i_Q 9g2 

NGC 4458 

EO-1 

16'= 

53 

0 

u.uuu/_o 0003 

0.0041°-°*° 

1.121 

- 

0 047+0-041 
U.U4/_0 047 

NGC 4459 

SAO 

7.1^ 

54 

1 

o-o7it°:°I^ 

n 1 81 +0.018 
U.i«i_o 017 

6.53 

0 c:7+U.U5 

U-O^_o 06 

0 477+0-023 
' -0.022 

NGC 4472 

E2 

7.9=^ 

51 

1 

0.085t°;«“ 

16 096'*'^'®^^ 

lO.UyO_0 330 

35.51 

1 nr+0.002 

-‘-•^9_0.002 

90 r^'7+0.592 

ZO.U/-0 007 

NGC 4473 

E5 

12.2™ 

51 

1 

0 074+0-038 

u.u/4_0 025 

0 0[5 c:+O.O 21 

U.UOO_o 017 

5.57 

0 7K+O.I9 

O'- 'O_0.20 

0 136+0-023 

U.iOU_o 024 

NGC 4477 

SBO 

9.6® 

55 

1 

0 030+0-021 
U.UOU_o 012 

0 740+0 027 

U. /4tu_0 027 

6.13 

0'-o'-L_o.i9 

1 (077+0.020 
‘ -0.020 

NGC 4494 

El-2 

12‘ 

0 

0 

7.413tJ;“ti 

0 014+0-012 
U.Ui^_0_oi4 

17.246 

0 7K+O.43 

O'- 'O_0.37 

u.zio_o 213 

NGC 4526 

SABO 

1.7*' 

18 

1 

0.0281°;“" 

0 506+0-109 

U.OUD_o_io9 

13.31 

^•^^-0.02 

9 9qo”I"0.234 

z.zoo_o 261 

NGC 4552 

EO-1 

9.6'= 

40 

1 

0 041+0 048 

U.U^l_0 022 

9 907+0-110 
z.zu/ _o 811 

7.83 

0 6+0-01 
0-O_o.oi 

1. <yU-o 368 

NGC 4564 

E 

5.9'= 

50 

1 

u.UiO_o_oo6 

0 00C+O.OO8 

U.UUO_o_oo5 

2.54 

- 

n 048+0-028 

U.U45_o 04g 

NGC 4596 

SBO 

11" 

36 

1 

0 0O7+O-OO5p 

u.uu^_o 003 

0 008+0-026 
U.Uy5_o 020 

5.67 

n 55+0-09 

0 399+0-050 

u.oyy-o 042 

NGC 4621 

El-2 

10.5j 

36 

1 

0 01Q+O-O32 

U.Uiy_o 012 

0.0531°;*°° 

5.54 

0 1 c+0-05 

y-lo_0.07 

0 362"*"^'^^^ 

U.90Z_o 

NGC 4636 

EO-1 

10.3 

21 

0 

u.uz^_0 009 

20.0281°;°°° 

10.84 

r» y^+0.002 

'4- <^_0.002 

-,K 19+0.190 

10.1Z_Q I^g2 

NGC 4649 

E5 

11.7" 

38 

1 

0 129+0 049 
u.izy_o 035 

11 090+^-^^^ 
ii.uyu_o 993 

29.68 

0 91+0-003 
^•^J-_o.o03 

y.zyZ-o_456 

NGC 4697 

E6 

8.9" 

0 

0 

0 072+0 042 

U.U/Z_o 027 

0 909+0.044 
u.zoz_o 036 

8.2 

0 39+0.04 

0 890+0-074 
U.52U_o Qgg 

NGC 4710 

SAO 

- 

8 

1 

u.ozo_o 070 

0.0871°-°;° 

4.45 

0 74+0.08 

U- /^_0.09 

u.z^y-0 023 

NGC 5422 

SO 

5* 

6 

0 

- 

0.017l°-°*4 

5.349 

n 15+0-80 

0 4tx8+0-163 

U.458_o_453 

NGC 5576 

E3 

10.2j 

6 

0 

0 01 3+0.008 
U.Ui-3_0_005 

0 036+0-042 

U.UOD_o 020 

7.672 

o'-^^-0.24 

0 1 qo+0.093 

U.iy»-o 066 

NGC 5813 

El-2 

18.3" 

23 

0 

u.uoz_o 019 

68.870i;-°“ 

22.443 

y+0.002 

'O- ‘ -0.002 

4R co+0.659 
4b.b5_Q og7 

NGC 5846 

EO-1 

13.5" 

18 

0 

n 087+0-045 

U.U«/_o_o30 

28.260i;-;™ 

23.8 

0 71+0.01 
'^•'-‘^-0.01 

20 37+0-417 
ZU.O/-0 424 

NGC 5866 

SAO 

1.8'* 

4 

0 

0 044+0-061 
U.U^4_0_o25 

0-4731°-°°; 

8.1 

rv 90 + O.O6 
^•^'3_o 03 

0 81 O+0.232 
U.5IO-0 193 

NGC 7457 

SAO 

2® 

- 

0 

u.uuz_o 001 

0.0031^.^3 

1.727 

- 

0 047+0-039 
U.U4/_0 047 


Note. — ★: morphology type taken from NED, not related to the parameter A/yT that we used to distinguish Es and SOs in this work. 
1 represents member galaxies of the Virgo Cluster; otherwise 0. f: derived within 2re. (a). Trager et al. (2000); (b). Crocker et al. 
(2012); (c). Terlevich & Duncan (2002); (d). Howell (2005); (e). Zhang et al. (2008); (f). Lees et al. (1991); (g). McDermid et al. (2006); 
(h). Gallagher et al. (2008); (i). Humphrey et al. (2006); (j). Idiart et al. (2007); (k). Sanchez-Blazquez et al. (2006); (1). Sil’Chenko (2006); 
(m). Kuntschner et al. (2010); (n). Sil’Chenko Chilingarian (2011); (o) Amblard et al. (2014); (p). Davis et al. (2014); (q). Shapiro et 
al. (2010). (*) Boroson et al. (2011). 
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malization of the apec thermal emission model, which is 
defined as 

where Da is the angular distance to the galaxy, and V is 
the volume of the source region. With the hot gas density 
and volume, we obtained the hot gas mass for each 

galaxy in the sample, as listed in Table 3. 

3.4. 2MASS photometry 

We used /f-band luminosities to characterize the stel¬ 
lar masses of the sample galaxies, since the iti-band ac¬ 
curately represents the largely old stellar populations of 
ETGs. These galaxies’ if-band luminosities were derived 
from Two Micron All Sky Survey {2MASS) (Skrutskie 
et al. 2006) archived images. The regions used for K- 
band photometry were the same as those used in the 
X-ray analyses. Bright nuclear and foreground sources 
(detected by eye) were excluded and refilled with a local 
surface brightness component using the dmfilth tool in 
CIAO 4.5. Source counts were obtained by subtracting 
the local background component. We converted source 
counts to the corresponding magnitude and corrected for 
Galactic extinction. The iiT-band solar luminosity was 
assumed to be Lkq = 5.67 x 10^^ ergs s“^ (Mannucci et 
al. 2005). The sample galaxies’ if-band luminosities Lk 
are listed in Table 2. 

3.5. Galaxy environment 

The Sloan Digital Sky Survey (SDSS) (York et al. 
2000) was used to assess the galaxy density around each 
galaxy in the sample. Using the main spectroscopic sam¬ 
ple from Data Release 7 (DR7, Abazajian et al. 2009), 
we count neighboring galaxies within a projected radius 
of 500 kpc and within a redshift difference of \z\ < 1500 
kms“^. Spectroscopic redshifts were used, since they 
are much more accurate than photometric redshifts. To 
compensate for the differences in the detection limit of 
galaxies at different distances, we included all galaxies 
brighter than = 17.8 for the most distant galaxies 
and used a correspondingly brighter magnitude thresh¬ 
old for more nearby galaxies. The number of galaxies in 
the neighborhood of each sample galaxy is listed in Ta¬ 
ble 3. Galaxies with the most numerous neighbors tend 
to reside in the Virgo Cluster. Four out of 42 galaxies 
in our sample are not well covered by the SDSS and are 
described in §2. 

We supplemented this environmental study of 42 galax¬ 
ies possessing sufficiently deep, high resolution Chan¬ 
dra X-ray data by investigating a much larger sample 
of galaxies observed with the (lower resolution) ROSAT 
X-ray telescope. The O’Sullivan et al. (2001) ROSAT 
galaxy survey includes 430 ETGs, 209 of which were de¬ 
tected. The histogram in Figure|^shows the distribution 
of the number of neighboring galaxies around each of the 
430 ETGs in the O’Sullivan et al. (2001) sample. The 
galaxy sample was taken from NEEIdand neighbors were 
defined to be those galaxies within a 500 kpc projected 
radius of each ROSAT galaxy. Despite the statistical 
unevenness of both the NED and ROSAT samples, there 

® http://ned.ipac.caltech.edu/ 



Fig. 2.— Histogram of the number of neighbor galaxies within 
500 kpc of each galaxy based on RC3 for 430 galaxies taken from 
O’Sullivan et al. (2001). We take 15 nearby neighbor galaxies as 
a natural cut between galaxies in dense environments and galaxies 
in relatively isolated environment. 


45 

44 

43 

42 


W) 

_o 

^ 40 
39 
38 
37 

8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 

logL^ 

Fig. 3.— The total X-ray luminosity Lx (including LMXBs and 
stellar sources) as a function of Lg for 209 ETGs in isolated and 
dense environments taken from the ROSAT survey of O’Sullivan 
et al. (2001). [see the electronic edition of the journal for a color 
version of this figure.] 

seems to be a natural division between high-density and 
low-density environments at the level of ~ 15 neighbors. 
Of the 209 ETGs detected in the O’Sullivan et al. (2001) 
ROSAT survey, 146 are in dense environments, while 63 
are in low-density environments. Figurej^plots the X-ray 
and optical (blue) luminosities of the 209 ROSAT detec¬ 
tions; the galaxies are labelled by environment, with 15 
neighbors dividing high-density from low-density envi¬ 
ronments. The high- and low-density samples are almost 
completely overlapping. 

4. RESULTS 

4.1. The scatter in the Lv -Lk relation 

-^gas 

The gaseous X-ray luminosities of our sample 

galaxies are plotted against the galaxies’ stellar AT-band 


• galaxies in the field 
A galaxies in groups and clusters 


• A. 















luminosities Lk in Figure]^ ranges from ~ 3x 10^^ 
to 7 X 10^^ ergs s“^, while Lx ranges from ^ 10^° 
to 3 X 10^^ Lxq (see Table 3). The scatter in the 
Lx^^-Lk correlation ranges up to a factor of 1000, 
demonstrating the diversity of galaxies in our sample. 
We fit the Lx^^^-Lk relation to a single power law, 
log(Lxg„) = A log(Lif) + B, and obtained best-fit val¬ 
ues for the slope ofA = 2.3±0.3 and the intercept of 
B = 14.6±3.3, as indicated by the solid black line in Fig¬ 
ure The slope we obtained is slightly shallower than 
(but within the errors of) that obtained for a Chandra 
sample of 30 ETGs by Boroson et al. (2011), who found 
a best-fit slope of 2.6 ± 0.4. Our sample contains rela¬ 
tively more faint galaxies, which may have a shallower 
Lx^^^-Lk slope than the sample as a whole. Moreover, 
our study uses the 0.1-2 keV energy band, rather than 
the 0.5-2 keV band used by Boroson et al. (2011). This 
would flatten our slope relative to the 0.5-2 keV slope, 
since X-ray fainter galaxies tend to have lower tempera¬ 
tures, thus have a higher fraction of their total gaseous 
emission below 0.5 keV than do more X-ray luminous 
galaxies. To test this, we re-performed our analysis us¬ 
ing the 0.5-2.0 keV band and obtained a slope of 2.4±0.3, 
in better agreement with the Boroson et al. (2011) result. 
Another major difference between our work and that of 
Boroson et al. (2011) is that we use metric radii (two op¬ 
tical effective radii) for the extraction regions in both 
the X-ray and AT-band analyses, while Boroson et al. 
(2011) use background-limited photometric radii (radii 
where the diffuse emission reaches the background level) 
in their X-ray analysis and their Lx values were taken 
from NED. This could also cause differences in our re¬ 
sults; nevertheless, our approach is more self-consistent. 
In the Appendix, we further examine the impact of vari¬ 
ations in the best-fit Lx^^^-Lx relation (including using 
a best-fit broken power law) upon our results. 

In Figure we show the Lx^^^-Lx relation for each 
of the three sets of sub-groupings (E/SO, mass, envi¬ 
ronment) described in §2. We find that massive galax¬ 
ies, “true elliptical” galaxies, and galaxies in groups and 
clusters generally belong to a single population (denoted 
population A), while low-mass galaxies, “true lenticular” 
galaxies, and field galaxies generally belong to a separate 
population (denoted population B). 

If ETGs retained all their accumulated stellar-mass 
loss, the hot gas mass Mx ^., should in principle be lin¬ 
early related to Lx- Figure^shows Mx^^^ as a function 
of Lx for our sample. The scatter of the Mx^^^-Lx re¬ 
lation is smaller than the Lx ^^-Lx relation (Figure |4). 
This is manifested in the smaller uncertainties in its bes - 
fit power fit: a best-fit slope of A = 1.84 ± 0.17 and in¬ 
tercept B = —12.1 ± 1.8, as indicated by the black solid 
line in Figure This result, although steeper than 

linear, slightly reduces the discrepancy we had between 
the hot gas content and the optical light of ETGs. We 
further compared the total gas mass Mgas (hot gas plus 
cold gas: Mgas = AIx^^^ + Mhi + Mua) as a function of 
Lx for our sample as shown in Figure fright. We ob¬ 
tained a best-fit slope of A = 1.11 ± 0.21 and intercept 
B = —3.61 ±2.20. The scatter of this Mgas —T/c relation 
is larger than in the relation between the hot gas mass 
and Lx- This suggests that the cold gas content is un¬ 
likely to be from accumulated stellar-mass loss or cooling 


from hot gas. However, we note that the total gas mass 
Mgas increases almost linearly with Lx, which is exactly 
what we expect if stellar-mass loss is the primary source 
of all gaseous components. 

In subsequent sections (§4.2-§4.8), we relate the 
Lx — Lk relation to various internal and external fac- 
tors: galaxy age, cold gas masses (A/hi, Mua, Mhi ± 
Mhj), star formation rate (SFR), total galaxy mass 
(Mtot, Te, Mtot/re), hot gas temperature {Tx), stellar ve¬ 
locity dispersion (u), a rotational parameter (A), galaxy 
eccentricity (e), their combination (A/-\/e), and the local 
density of galaxies in the neighborhood of each galaxy. 
We also investigate the deviations in Txg„ from the best- 
fit Lx — Lx relation and assess whether these devia- 
tions [defined as ALx^^, = log(Lxg„) - 2.31og(Lx) - 
14.6] are correlated with the internal and external fac¬ 
tors listed above. We use the Spearman correlation in the 
Penn State statistical package ASUKVnto quantify the 
correlations of these relations, which allows us to prop¬ 
erly treat any upper limits. The correlation coefficient 
(with uncertainties) and its associated null hypothesis 
probability for any relation between Lx^^^/Lx and each 
factor are listed in Table 4. The results for the relation 
between the residuals ALxgaa and each factor are listed 
in Table 5. We also repeated all these analyses using 
Mxg„^ jLx and analogous deviations AMxg,,,, instead of 
Lx^^^/Lx and ALxgasi obtained very similar results, 
so they are not presented here. 

4.2. Stellar age 

Lxg^jLx and ALx^^^ are plotted against the st ellar 
ages of our sample galaxies in Figures Al-top-left and |A4[ 
respectively. These figures show that galaxies with older 
stellar populations tend to have relatively larger X-ray 
halos. Older galaxies have a longer time to accumulate 
hot gaseous halos. This trend was suggested by Boroson 
et al. (2011) for nearby ETGs. In contrast, Givano et 
al. (2014) found in the Chandra GOSMOS survey that 
younger galaxies tend to have larger Lx^^^/Lx- This 
may be related to younger galaxies experiencing recent 
major mergers. Note that galaxies in the Chandra COS¬ 
MOS survey extend out to z = 1.5, so may have different 
evolutionary histories compared with nearby ETGs. At 
the same time, the galaxy ages listed in Table 3 are taken 
from heterogenous sources. The lack of consistency be¬ 
tween these analyses and the fact that many galaxies 
show radial age gradients may further complicate this 
study. We also note that the correlation between age and 
ALxga,= seems to only exist among population A galaxies. 

We estimated the stellar-mass loss rate for each galaxy 
by using the current value determined for elliptical galax¬ 
ies by Knapp et al. (1992): ^0.0021 Lx/Lxq Mq Gyr“^. 
This can be regarded as a lower limit to the average 
stellar-mass loss rate over cosmic time. We calculated 
the expected gas mass for each galaxy by multiplying 
its stellar-mass loss rate by its stellar age. In Figure 0 
we compare the observed total gas mass (hot gas plus 
cold gas) of each galaxy to the expected gas mass (the 
solid black line indicates equality). For most galaxies, the 
observed gas mass is much less than that expected from 
accumulated stellar-mass loss. This discrepancy suggests 

® http://www2.astro.psu.edu/statcodes/asurv 



9 


OX) 

S 


OX) 

o 



10.0 10.5 11.0 11.5 


0C1O23^N4261 
4^4821 ★N4278 
AN1023 ♦N4342 

★n1266 HN4365 

♦N2768 #N4374 
BM2778 4iN4382 
®N3377 ^N4406 
4^^3379 ★N4458 
AN3384 ♦N4459 
★N3599 BN4472 
♦N3607 #N4473 
B43608 4’N4477 
#43665 ^N4494 
4<^4203 ★N4526 


♦ N4552 
■N4564 
#N4596 
4'N4621 
AN4636 
★N4649 

♦ N4697 
■N4710 
#N5422 
4'N5576 
AN5813 
★N5846 

♦ N5866 
HN7457 


12.0 


12.5 


Fig. 4.— Gaseous X-ray luminosities as a function of X-band luminosities for all galaxies in this work. The black solid line is the best-fit 
Lx,^^-Lk relation ( log(LxgaJ = 2.31og(LK) + 14.6). 


that gas has been driven out of these galaxies. Alterna¬ 
tively, this may indicate that gas has been distributed to 
larger radii than our extraction apertures and is more ex¬ 
tended than the stellar light distributions. Humphrey et 
al. (2011) studied the isolated elliptical galaxy NGC 720 
in great detail using deep Chandra and Suzaku observa¬ 
tions, the combination of which provides both high spa¬ 
tial resolution and low instrumental background. They 
detected hot ISM emission to beyond 50 kpc and extrap¬ 
olate its hot gas properties out to the virial radius 
300 kpc). Based on its measured and extrapolated 
gas mass profiles, its enclosed gas mass increases by 10 
times (600 times) out to 45 kpc (300 kpc) compared to 
its gas mass within 2 re (~10 kpc). If we scale the gas 
mass of galaxies in our sample out to 60 kpc and 300 kpc, 
based on the gas mass profile of NGC 720, the equality 
between the observed total mass and the expected gas 
mass would shift to the red and blue dashed lines, re¬ 
spectively, in Figure In principle, we may be seeing 
the inner parts of the history of accumulated stellar-mass 
loss, while the bulk has been pushed beyond the optical 
extent of the galaxies. 

4.3. Cold gas content and star formation 

More than 20% of the ETGs in the ATLAS^^ multi¬ 
wavelength survey are observed to contain atomic and 
molecular gas in significant amounts (at least 10 ^- 
10® Mq). Su & Irwin (2013) proposed that the hot ISM 
of ETGs may have accreted cold gas that was subse¬ 
quently heated to an X-ray emitting phase. We thus ex¬ 
pect that the hot gas and cold gas content of ETGs may 
be related. We investigated the relation between the hot 
gas content in our sample galaxies and their cold gas 
masses: Mhi, and Mhi -f Mn^. We found that the 


deviations ALxgas generally increase with t he co ld gas 
masses in these galaxies, as shown in Figures ^ to |A 6 [ 
in particular, the hot gas in field galaxies is correlated 
with their molecular gas masses, as shown in Figure 
(p = 0.714 ± 0.117), while there is almost no correla¬ 
tion for population A galaxies. Figure^ top-right shows 
that H 2 is detected only in “true lenticular” galaxies in 
our sample. The effect of cold gas on galaxies in groups 
and clusters may be overshadowed by more complicated 
interactions with the ICM. 

These appreciable amounts of cold gas in ETGa are 
found to be associated with low levels of star formation, 
with rates ranging up to ~ lM 0 /yr (e.g. Temi et al. 
2009; Amblard et al. 2014). Diffuse hot gas can be pro¬ 
duced by the energy output associated with star forma¬ 
tion in late-type galaxies (e.g. Mineo et al. 2014). Stin¬ 
son et al. (2013) show that hot gaseous halos form at the 
same time as star formation develops. Current star for¬ 
mation in ETGs, even at low levels, could have affected 
the ISM hot gas content and contribute to the scatter in 
the Lxg^^-Lx relation. We found that is posi¬ 

tively correlated with the SFR, as shown in Figure A 
single obvious outlier is NGC 4494, an X-ray faint galaxy 
with a very large SFR of ~ 7.0 Mq/yt. NGC 4494 was 
also noted to have an unusually high SFR in our source 
reference, Amblard et al. (2014), who studied the cor¬ 
relation between SFR and dust content; Amblard et al. 
(2014) suggested that its anomalously high SFR may be 
related to its AGN and their SED fitting may have led 
to an unphysical solution. In fact, Wu & Gao (2006) and 
Satyapal et al. (2005) report a 50 times smaller SFR for 
NGC 4494, which is more consistent with its very small 
molecular hydrogen mass Mhj- Excluding NGC 4494, 
we found a very strong correlation between ALx and 
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Fig. 5.— Gaseous X-ray luminosities as a function of X-band 
luminosities for all galaxies in this work. The black solid line is 
the best-fit relation ( log(Z;Xgas) = 2.31og(Lx) + 14.6) 

for all sample galaxies. Dashed lines indicate the best-fit of each 
sub-group, top: “true lenticular” galaxies (brown) and “true ellip¬ 
tical” galaxies (green), middle: high mass galaxies (red) and low 
mass galaxies (black), bottom: field galaxies (blue) and galaxies in 
groups and clusters (magenta). 

SFR. 


4.4. Gravitational potential 

Galaxy mass has been regarded as one of the most 
plausible factors in influencing the hot gas content of 
ETGs (e.g. Mathews et al. 2006). Of all the correlations 


described above, we find that /Lk is most strongly 
cor relat ed with Mtot, the total mass within 1 (see Fig¬ 
ure A2-bottom left). Since Mtot is measured within Ire, 
it is dominated by the stellar mass, which scales with 
Lx (so does re). The measurement of Mtot is not within 
large enough radii to be dominated by dark matter. The 
Lxg^^-Lx relation is very steep, with Lx^^^ depending on 
Lk to a power of at least 2. This means the Lx^^^/Lk 
ratio has a residual dependence on Lk- Thus, the appar¬ 
ent Lx^^,/Lx-Mtot correlation could be a by-product of 
the intrinsic dependence of Lx^^^/Lx on Lx- 
We find that the deviations ALx^^^ also generally in¬ 
crease with Mtot (as well as re and Mtot/ re) fo r our sam¬ 
ple galaxies, as demonstrated in Figure |A7[ This may 
reflect that it is easier for more massive galaxies to retain 
their hot gas. A closer look shows that ALx^^^ increases 
with these factors only in more massive galaxies, while 
the deviations ALx^^^ are actually anti-correlated with 
Mtot in low mass galaxies. We show in the Appendix 
that this latter trend does not survive tests of system¬ 
atic errors. 


4.5. Hot gas temperature and stellar velocity dispersion 

The stellar velocity dispersion and the temperature of a 
hydrostatic gaseous halo in an ETG are partial measures 
of the depth of its gravitational potential in a thermalized 
system (the gas temperature can also be a sensitive in¬ 
dicator of non-gravitational energy feedback). Figure 10 
indicates that those population A galaxies with hotter 
atmospheres generally have larger ALx^.^^ - It is difficult 
to obtain accurate hot ISM temperatures in low mass 
galaxie^ Previous studies show that X-ray luminosities 
increase with ISM temperature in ETGs (Boroson et al. 
2011; Matthew et al. 2006; and Kim & Fabbiano 2013). 
Our study focuses on the hot gas content per unit stellar 
light, rather than the X-ray luminosity itself. We find 
that ALxg„ and Tx are positively correlated. The rela¬ 
tion between ALx^^^ and a is similar to the ALx^^^-Mtot 
relation. 


4.6. Rotation and ellipticity 

The ATLAS^^ project found that a parameter Xj\fe, 
combining a galaxy rotational parameter A with a galac¬ 
tic eccentricity e, distinguishes slowly-rotating “true el¬ 
liptical” from fast-rotating “true lenticular” galaxies 
(Emsellem et al. 2011) . The parameter A is robust mea¬ 
sure of a galaxy’s angular momentum and is defined as 

(Jesseit et al. 2009), where i?, V, and cr are the circular 
radius (1 Fe), circular velocity and stellar velocity disper- 
sion. The values of A and e of galaxies in our sample are 
taken from the ATLAS^^ survey (Emsellem et al. 2011) 
and listed in Table 2. We show as a function of 

^ The temperatures of NGC 821, NGC 3377, and NGC 3379 
were taken from Boroson et al. (2011), while we could not have 
their temperatures constrained in our fitting. This may be because 
Boroson et al. (2011) use a smaller Lx/^k for CV+AB component 
which leaves more counts for the hot gas components for these 
extremely faint galaxies. 
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Fig. 6.— left: hot gas mass as a function of /C-band luminosities for all galaxies in this work, right: total gas mass (the sum of hot gas 
and cold gas) as a function of i^-band luminosities for all galaxies in this work. The scatter of the total mass — Lk relation is larger than 
that of hot gas mass — Lk relation. Color code is the same as in Figure]^ 



Fig. 7.— Observed gas mass (the sum of cold and hot gas mass) 
as a function of expected gas mass (stellar age times stellar-mass 
loss rate). We use the average age (8.15 Gyr) for galaxies that do 
not have stellar age available in the literature and are marked by 
open symbols. Color code is the same as in Figure The solid 
black line indicates the equality between these two masses. For 
most galaxies, their observed gas mass is much less than expected. 
Dashed red line: equality between these two masses assuming gas 
distributes out to lOrg (~ 45 kpc). Dashed blue line: equality 
between these two masses assuming gas distributes out to rvir (~ 
300 kpc). [see the electronic edition of the journal for a color version 
of this figure.] 

the parameter X/^/e in Figure [h] We find that true el¬ 
liptical galaxies {X/y/e < 0.31) have a larger X-ray excess 
than true lenticular galaxies (X/y/e > 0.31). 


We also compared the hot gas content to the rotation 
parameter A and ellipticity e separately. We find that the 
deviations ALy are anti-correlated with rotation and 
flatness, as shown in Figures |A8| and |A9[ respectively; 
however, this only holds for population A galaxies. The 
situation is further complicated by the fact that flatness 
and rotation are usually highly correlated, in the sense 
that only flat galaxies can support fast rotation. This 
is the case for galaxies in our sample. It is difficult to 
tell whether the above two correlations have real physi¬ 
cal origins or one correlation is simply a byproduct of the 
other. The best way to test this is to study the correla¬ 
tion of hot gas and rotation for galaxies that have similar 
shapes or to study the correlation between hot gas and 
ellipticity for galaxies that have similar rotation. 

Based on the data in Table 2 and Figures |X8] and [XO} 
we find that A = 0.15 provides a natural cut separating 
slow-rotating and fast-rotating galaxies and e = 0.35 pro¬ 
vides a natural separation between flat and round galax¬ 
ies. We divide our sample into flat (e > 0.35) and round 
(e < 0.35) galaxies in the ALxgaE>“™t^tion relation, as 
shown in Figure \T^left. Neither sub-group exhibits a 
significant correlation. As shown in Figure \12\ right, we 
divide our sample into fast-rotating (A > 0.15) and slow- 
rotating (A < 0.15) galaxies in the -ellipticity 

relation. Slow-rotating galaxies show a significant cor¬ 
relation between the negative residuals of ALxgas 
ellipticity. This result indicates the relatively important 
role of ellipticity. 

4.7. Environmental density 

It is debated whether intracluster gas has a net posi¬ 
tive or negative effect on the X-ray luminosities of cluster 
galaxies (Sun et al. 2007; Mulchaey & Jeltema 2010). For 
galaxies in dense environments, gaseous outflows can be 
suppressed by ICM pressure confinement. In more iso¬ 
lated environments, gaseous outflows could be respon- 
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Fig. 8.— ALxgas ^ function of molecular gas mass, top-left-, color code is the same as in Figurej^ top-right “true lenticular” galaxies 
(brown) and “true elliptical” galaxies (green), bottom-left high mass galaxies (red) and low mass galaxies (black), bottom-right: field 
galaxies (blue) and galaxies in groups and clusters (magenta). 


sible for the scatter in Lx^^^/Lk for ETGs. From the 
SDSS we obtained the number of neighboring galaxies for 
38 out of 42 galaxies in our sample, as listed in Table 3. 
Figure AlO shows Lx^^jLx and ALx^^^ as a function of 
the number of nearby galaxies. Galaxies in higher galaxy 
density environments contain more hot gas per unit stel¬ 
lar light (although this trend seems very mild), as found 
by Brown & Bregman (2000) using ROSAT and ASCA. 
This may be a result of pressure confinement of the ICM, 
or the accretion of the surrounding IGM into ISM (Pinino 
et al. 2005), or they may be massive slow-rotators, 


4.8. Ram pressure stripping in the Virgo Cluster 


The Virgo Cluster is the nearest relaxed galaxy cluster, 
residing at a distance of only 16 Mpc. Gen et al. (2014) 
show through simulations that ram pressure stripping of 
galaxies starts to become effective within 3 Mpc of the 
cluster center. In our sample, 15 galaxies have their de- 
projected radii within 3 Mpc of M87, the central galaxy 
of the Virgo Cluster. This allows us to study the effects 
of ram pressure on many galaxies in a single cluster. A 
quantitative description of the strength of ram pressure 
stripping was first presented by Gunn & Gott (1972). 
Based on a static force argument, a stripping condition 
can be derived for a gaseous disc moving face-on through 
the ICM. Hot gas in the ISM will be stripped when the 
ram pressure (Pram = Picm'*^^) exceeds the gravitational 
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SFR (MJyr) SFR (MJyr) 

Fig. 9.— ALjfgas 3,s a function of star formation rate, top-left: color code is the same as in Figurej^ top-right: “true lenticular” galaxies 
(brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black), bottom-right: field 
galaxies (blue) and galaxies in groups and clusters (magenta). 


restoring force per unit area, Pgrav- McCarthy et al. 
(2008) developed an analogous model for the ram pres¬ 
sure stripping of galaxies with a spherically-symmetric 
gas distribution. Their model, which is more suitable for 
ETGs with an extended hot gas halo rather than a cold 
gaseous disk, yields 

D F GifftotPiSM /T'l 

.f^ram ^ o p > 

where Mtot is the total mass of the galaxy, is the 
galaxy ISM gas density, and is the radius of the 
galaxy at which the stripping occurs. We consequently 
obtained the ratio of ram pressure to gravitational restor¬ 
ing pressure: 


Rrara _ ^ PiCM ^ -^ISM ^2^ 

Cgrav PiSM Gilftot 

We adopt eq.(2) to calculate the instantaneous 
strength of ram pressure stripping relative to the gravita¬ 
tional restoring force. We estimate Pj^M from the depro- 
jected gas density profile of the Virgo Cluster in the ra¬ 
dial range of 0.3-1 Tvir, as derived from 13 XMM-Newton 
observations out to r^ir (1.2 Mpc). This results in a de- 
projected density profile of the form pj(g„ oc r~^-^ (Ur¬ 
ban et al. 2011), where r is the 3D separation between a 
member galaxy and M87, as taken from ATLAS^'^. The 
relative velocity v is chosen to be VS times the observed 
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Fig. 10.— ALxgaa 3-® ^ function of hot ISM temperature, top-left: color code is the same as in Figure top-right: “true lenticular” 
galaxies (brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black), bottom-right: 
field galaxies (blue) and galaxies in groups and clusters (magenta). 


radial velocity relative to M87, obtained from NED, and 
is chosen to be 2re. We obtained values of 
through the best-fit norm s of the apec component in our 
spectral analysis. Figure 13 shows the ratio Lx,^^/Lk 


and ALx^^ as a function oTthe ratio Pram/^grav for the 
15 Virgo Cluster ETGs in our sample, as described in §2. 
X-ray faint galaxies tend to have larger Pram/J^grav than 
X-ray luminous galaxies. This trend implies that X-ray 
luminous galaxies have retained a larger fraction of their 
gas against ram pressure, while X-ray faint galaxies have 
lost a larger fraction of their gas to ram pressure. 

5. DISCUSSION 


The large scatter (up to a factor of 1000) in the 
Lk relation for ETGs has been a long-standing puzzle 


in extragalactic astronomy. We use Lx^^JLk and the 
residuals ALx^.^^ to describe the hot gas content per unit 
stellar light in our sample of 42 galaxies and compare 
them to various internal and external factors. We find 
that the X-ray luminosity deviations ALx^^^ are most 
strongly correlated with the galaxies’ star formation rates 
and hot gas temperatures. We also find two distinct pop¬ 
ulations of ETGs. These results point towards some im¬ 
portant implications that are discussed in the following 
subsections. 

5.1. Population A and population B galaxies 

Galaxies in our sample tend to belong to two different 
populations, as shown in Figure Population A galax- 
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TABLE 4 

Summary of correlation of with each factor 


Relations 

Correlation coefficient p 

Errors (STA) 

Null hypothesis probability 

Number of galaxies 

Correlation coefficient r 

Age 

0.439 

±0.034 

0.0084 

37 

0.215 

M(H i) 

0.412 

±0.033 

0.0100 

39 

0.182 

M(H2) 

0.152 

±0.036 

0.3310 

42 

0.098 

M(H2) -l-M(H i) 

0.336 

±0.030 

0.0314 

42 

0.224 

SFR* 

0.568 

±0.068 

0.0006 

38 

0.425 

Tx 

0.515 

±0.184 

0.0023 

36 

0.352 

(7 

0.553 

±0.025 

0.0004 

42 

0.317 

Mtot (r < r-e) 

0.639 

±0.070 

0.0000 

42 

0.457 

re 

0.622 

±0.024 

0.0001 

42 

0.460 

^tot/re 

0.611 

±0.027 

0.0001 

42 

0.171 


-0.441 

±0.026 

0.0047 

42 

-0.285 

Rotation A 

-0.439 

±0.028 

0.0050 

42 

-0.282 

Ellipticity e 

-0.280 

±0.034 

0.0732 

42 

-0.193 

SDSS 

0.346 

±0.033 

0.0352 

38 

0.185 

-Eram/-^gra 

-0.479 

±0.041 

0.0733 

15 

-0.333 


Note. — * not including NGC 4494. 


TABLE 5 

Summary of correlation oe ALxg ^^ with each factor 


Relations 

Coefficient p'f’ 

Errors (STA, SYS#) 

Null hypothesis prob. 

No. of galaxies 

Correlation coefficient 

Age 

0.283 

±0.080, 

o 

o 

d> 

1 

±0.009, 

-0.041 

0.0892 

37 

0.182 

M(H i) 

0.319 

±0.067, 

±0.049, 

-0.067, 

-0.028 

0.0495 

39 

0.153 

M(H2) 

0.213 

±0.058, 

±0.020, 

1 

O 

b 

to 

±0.008 

0.1717 

42 

0.128 

M(H2) ±M(H i) 

0.296 

±0.060, 

±0.032, 

-0.041, 

±0.057 

0.0584 

42 

0.227 

SFR* 

0.536 

±0.103, 

±0.025, 

-0.065, 

-0.026 

0.0011 

38 

0.428 

Tx 

0.476 

±0.186, 

1 

o 

b 

to 

±0.004, 

-0.098 

0.0049 

36 

0.337 

<J 

0.220 

±0.063, 

±0.106, 

-0.131, 

-0.088 

0.1595 

42 

0.064 

Mtot {r < To) 

0.175 

±0.065, 

±0.127, 

-0.149, 

-0.062 

0.2631 

42 

0.120 

re 

0.218 

±0.064, 

±0.114, 

-0.141, 

-0.070 

0.1631 

42 

0.151 

A^tot/^e 

0.231 

±0.070, 

±0.117, 

-0.139, 

-0.076 

0.1396 

42 

0.171 

A/v/5 

-0.217 

±0.070, 

-0.060, 

±0.061, 

±0.101 

0.1656 

42 

-0.138 

Rotation A 

-0.247 

±0.071, 

-0.064, 

±0.064, 

±0.106 

0.1140 

42 

-0.159 

Ellipticity e 

-0.272 

±0.063, 

±0.009, 

±0.011, 

±0.060 

0.0821 

42 

-0.261 

SDSS 

0.362 

±0.063, 

-0.011, 

±0.006, 

-0.086 

0.0276 

38 

0.202 

-Eram/-Egra 

-0.343 

±0.073, 

-0.032, 

±0.100, 

±0.168 

0.1995 

15 

-0.295 


Note. — 

respectively in Figure 
Kendall’s tau correlatidnT 


not i nclud ing NGC 4494. Systematic errors obtained using dashed red, blue, and black (broken power law) lines 
All as the benchmark of the Lx„^^lI jk relation, t determined by Spearman’s correlation, t determined by 


ies are generally X-ray luminous, massive, slowly rotat¬ 
ing, hot, big, and tend to reside in high galaxy density 
environments. Population B galaxies tend to be X-ray 
sub-luminous, low mass, fast-rotating, cool, small, and 
reside in the low galaxy density environments (the field). 
These factors are highly degenerate and it is difficult to 
pin down the driving factor. Perhaps galaxies in groups 
and clusters have experienced more mergers, which could 
slow down their rotation and produce more massive and 
bigger galaxies with deeper gravitational potentials. 

The best-fit power law characterizing each sub-group is 
listed in Table 6. The slopes of the sub-group Lx^^^-Lk 
relations are somewhat steeper for population A galax¬ 
ies than for population B galaxies; however, the slopes 
are significantly different for only the high-mass / low- 
mass sub-groups (which differ by more than 3 ct). We 
define the dispersion in — Lk for each sub-group 

as \J ^ ALxgaa^j where ALx^^^ is the difference be¬ 


tween the observed Lx and the best-fit Lx -Lx rela- 
tion for each sub-group and N is the number of galaxies 
in each sub-group, as listed in Table 6. Galaxies in the 
field have significantly larger dispersion than galaxies in 
groups and clusters (by 20%). This is consistent with our 
investigation of the environments of 209 ETGs detected 
by ROSAT, as shown in Figure 
Galaxies in groups and clusters experience complicated 
environmental effects. Relatively lower mass galaxies 
with a lower ISM temperature (thus a small ISM pres¬ 
sure) may suffer from ram pressure stripping and their 
atmospheres may evaporate due to the heating of the sur¬ 
rounding ICM. The atmospheres of massive galaxies are 
relatively hot with a high gas pressure. The surrounding 
ICM may provide pressure confinement. These processes 
make faint galaxies fainter and luminous galaxies more 
luminous, leading to a steeper slope in the Lx^^^-Lk re¬ 
lation for Population A galaxies. 
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TABLE 6 

Best-fits of the Ly -Lr- to a single power, law for each subgroup and their dispersions 

gas 


Subgroup 

Slope 

Difference in slopes 

Dispersion 

Number of galaxies 

Slope-Kelly* 

True elliptical 

2.545±0.508 


0.646 

12 

2.459 ± 0.772 

True lenticular 

1.619±0.373 


0.717 

30 

1.591 ± 0.342 

High mass 

3.222±0.700 

3 1 (T 

0.690 

20 

2.913 ±0.681 

Low mass 

0.843±0.394 


0.685 

22 

0.690 ±0.611 

Goups Sz Clusters 

2.515±0.326 

2 4 (T 

0.589 

21 

2.393 ±0.412 

Field 

1.374±0.481 


0.762 

21 

1.422 ±0.494 


Note. — * Slopes determined by the maximum likelihood procedure of Kelly (2007). 



Fig. 11.— ALxgaB function of the parameter A/-yfi. Color 
code is the same as in Figure The red dashed line indicates 
A/y^ = 0.31. We call galaxies with A/yfe < 0.31 true elliptical 
galaxies and galaxies with A/y^ > 0.31 true lenticular galaxies. 

5.2. Cold gas and stellar feedback 
5.2.1. Feedback from mild star formation 

In our study, we find that the hot gas content is most 
strongly correlated with SFR compared to all other fac¬ 
tors we have considered. This applies to both population 
A galaxies and population B galaxies. There are three ex¬ 
planations for this correlation. One explanation is that 
the hot gas transforms into cold gas through radiative 
cooling and forms into stars. We have three lines of evi¬ 
dence to rule out this scenario: ( 1 ). we hnd that the ra¬ 
diative cooling rate is smaller than the SFR for more than 
half of the galaxies. Thus, radiative cooling is unable to 
sustain the star formation. ( 2 ). ALx is positively cor¬ 
related with hot gas temperature and the SFR is mildly 
positively correlated with the hot gas entropy. Hot gas 
is not cooling fast. (3). the correlation between 
and SFR is stronger than that between ALx^^^ and Mhj . 
The second explanation is that stellar feedback energizes 
the ISM. Given that the star formation rates in our sam¬ 
ple galaxies are very low (mostly < O.lM 0 /yr) it seems 
unlikely that such star formation could drive significant 
amounts of gas out of the host galaxies; associated super¬ 


nova explosions (1 SN per 100 years per Mq/ji of star 
formation) would just heat up more gas into an X-ray 
emitting phase. Energetic feedback is expected to eas¬ 
ily sweep out the gaseous component in dwarf galaxies 
(Hopkins et al. 2012). ETGs in our sample are at least 
as massive as the Milky Way. Their gravitational poten¬ 
tials may be deep enough to retain their hot atmospheres. 
Thus, our result provides observational support of star 
formation feedback in early-type galaxies. The third ex¬ 
planation is that both low level star formation and en¬ 
hanced X-ray luminosities may both be consequences of 
recent wet minor mergers, as suggested by Civano et al. 
(2014). These processes need to be studied in more detail 
(especially via simulations) in the future. 

5.2.2. Origins of cold gas 

In this work we find that the hot gas content of pop¬ 
ulation B galaxies is significantly correlated with their 
neutral gas, particularly molecular gas mass. If the ra¬ 
diative cooling of hot gas is the main supplier of the 
cold gas, this would naturally explain why the cold gas 
content is observed to increase with the hot gas con¬ 
tent. In order to probe this scenario, we calculated the 
classical mass deposition rate for galaxies in our sam¬ 
ple: Me = 2gLuipLx/5kT (David et al. 2014). Then, 
we obtained the time required for each galaxy to form 
their cold gas contents through this deposition. In Fig¬ 
ure [H we compared the time needed to form the ob¬ 
serve amount of molecular gas and the stellar age of 
each galaxy. There is not enough time for field galaxies 
in our sample to have their cold gas supplied by hot gas 
cooling. Their cold gas may be acquired externally from 
mergers/accretions. However, it is possible for galaxies 
in groups and clusters to have most of their cold gas sup¬ 
plied by hot gas cooling. Sun et al. (2007) observed a con¬ 
nection between the hot coronae and the radio emission 
in cluster galaxies, which suggests that radiative cooling 
of the coronal gas may fuel the central black holes in 
galaxies in a cluster environment where galactic cold gas 
is difficult to acquire. 

Su & Irwin (2013) demonstrated an anti-correlation 
between molecular gas mass and the hot gas metallic- 
ity. One explanation for this trend would be that hot 
ISM is diluted by mixing with less metal-enriched molec¬ 
ular gas. This process would also add more hot gas to 
the ISM and lead to the ALx^^^-Mh^ correlation ob¬ 
served in this study. This is only effective in population 
B galaxies. Small mass systems have a smaller hot gas 
halo and they are easier to dilute. Galaxies in the field 
may have acquired cold gas through external processes 
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Fig. 12.— left ^ function of rotation A for flat galaxies (red) and round galaxies (blue), right ALxg^s ^ function of 

ellipticity e for fast-rotating galaxies (gray) and slow-rotating galaxies (yellow). 
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Fig. 13.— Lx^^^ILk and ALxg^a galaxies in the Virgo Cluster in this work as a function of the ratio of ram pressure to gravitational 
restoring pressure. Color code is the same as in Figure]^ [see the electronic edition of the journal for a color version of this figure.] 


such as recent mergers; these mergers would increase hot 
gas content. In contrast, galaxies in groups and clusters 
are more likely to have molecular gas supplied by stellar- 
mass loss or condensation from hot gas and their hot gas 
content may also be influenced by environmental factors 
such as ram pressure stripping. 

5.2.3. Star formation efficiency 

Star formation is observed to qualitatively correlate 
with its fuel, the cold gas. However, the exact process 
of star formation is still unclear. Various factors, such 


as environment, depth of potential well, metallicity, and 
AGN feedback, may have played noticeable roles in reg¬ 
ulating the star formation. The star formation efficiency 
(the ratio of SFR to cold gas mass) can vary dramatically 
among galaxies. In Figure]^ we show the relation be¬ 
tween SFR and cold gas mass for galaxies in our sample. 
SFR increases with cold gas mass for population B galax¬ 
ies, while population A galaxies tend to have constant 
SFR regardless of their cold gas mass. Various mecha¬ 
nisms may have interfered with the star formation pro¬ 
cess in population A galaxies. On the positive side, the 
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Fig. 14. — Time needed to form the observed amount of molecu¬ 
lar gas from radiative cooling of hot gas as a function of stellar age 
for field galaxies (blue) and galaxies in groups and clusters (ma¬ 
genta). The solid red line indicates the equality between these two 
time scales. We use the average age (8.15 Gyr) for galaxies that 
do not have stellar age available in the literature and are marked 
by open symbols. We use the average temperature (0.56 keV) for 
galaxies that do not have a constrained ISM temperature measure¬ 
ment and are marked by open symbols, [see the electronic edition 
of the journal for a color version of this figure.] 

ICM confinement may increase the gas density of cluster 
galaxies, which would boost star formation. Meanwhile, 
cluster galaxies may have their cold gas removed quickly 
through ram pressure stripping or heat conduction and 
appear to have a large SFR per unit gas. Moreover, 
star formation efficiency is found to be positively cor¬ 
related with gas metallicity (Yates & Kauffmann 2014). 
We demonstrated that cold gas in cluster galaxies and 
field galaxies may have different origins. Cold gas in clus¬ 
ter galaxies, formed through hot gas cooling, is probably 
contaminated, while cold gas in held galaxies, acquired 
externally, is more likely to be pristine. This could en¬ 
hance the star formation efficiency in cluster galaxies. On 
the negative side, Davis et al. (2014) demonstrated that 
the depletion time of cold gas positively correlates with 
the mean shear rate (A/0; A is the hrst Oort constant 
and O is the angular velocity). “True elliptical” galaxies 
rotate slower (high shear rates) and thereby may take 
longer time to form stars. Overall, the expected SFR 
and cold gas mass relation may have been washed out 
for population A galaxies due to various mechanisms. 

5.3. Galaxy masses 

Galaxy masses, dominated by dark matter halos, have 
been regarded as the most crucial factor in regulating 
the hot gas content of ETGs. Galactic winds supply 
gas to the hot gas halo; the explosion of supernovae 
provides energy. These processes would make massive 
galaxies brighter by adding more hot gas. However, 
these same processes may instead make low mass galax¬ 
ies fainter if the ejected gas becomes energized enough 
to escape the galaxy due to its shallower gravitational 


potential well. Not only strongly motivated in theory, 
the role of dark matter halos in retaining hot gas has 
been supported by observations of ETGs. Kim & Fab- 
biano (2013) demonstrated that Lx is highly corre¬ 
lated with the total mass for a small sample of ETGs, 
although the absolute values of Lx^^^ is not our primary 
interest. More relevantly, Mathews et al. (2006) found 
that Lx /Lk and even ALx^^^ by their definition (which 
is similar to ALx^^^ defined in our work) are correlated 
with their total masses. Galaxies in their sample are 
mostly bright elliptical galaxies with Lk > 2 x 10^^ Lkq 
and Lx > 3 x 10"^° ergs“^, with most residing at group 
centers. 

We do not observe a significant relation bet ween 
ALxg^^ and Mtot for galaxies in our sample (Eigure A7). 
The total mass used in Mathews et al. (2006) is the N-ray 
hydrostatic mass derived under the assumption of hydro¬ 
static equilibrium and spherical symmetry, extending out 
to fairly large radii. The total mass used by Kim & Fab- 
biano (2013) is taken from Deason et al. (2012) which is 
the dynamical mass within 5 . Both of these are very 

good proxies for dark matter mass. Unfortunately, the X- 
ray technique requires a large number of counts and can 
be very observationally expensive, in particular for rela¬ 
tively faint galaxies. The total mass estimated through 
stellar dynamics, extending to large radii, is not readily 
available since accurate mass estimates from stellar dy¬ 
namics are generally restricted to the very central regions 
of galaxies. The total mass used in this study is measured 
within 1 Te by the ATLAS^° survey, so this total mass is 
mostly made up of stellar mass and can not fairly reflect 
dark matter content. In order to extend the total mass 
within Ire out to 5re, we parameterize the dynamical 
mass of 11 ATLAS^^ galaxies studied by Deason et al. 
(2012) out to 5 re as a function of their total mass within 
1 re, as estimated in the ATLAS^^ survey. Applying this 
parameterization to all galaxies in our sample we found 
that the correlation between this scaled mass and ALx 

-^gas 

is still not significant. This does not rule out the role of 
a dark matter halo in regulating the hot gas content in 
ETGs. In some sense, hot gas temperature may serve as 
a proxy of potential well and we did observe that ALx^^^ 
increases with the ISM temperature in spite of the larger 
scatter. 

The lack of significant correlation may also be related 
to a complicated feedback process, in particular related 
to AGN activities. It has been well established that the 
bulge mass of a galaxy is highly correlated with the mass 
of the central supermassive black hole (e.g. Bentz et al. 
2008). Thus, more massive galaxies may have more ac¬ 
tive and stronger AGN outbursts, leading to stronger 
gaseous outflow. As a result, the gas content of more 
massive galaxies cannot be simply correlated with the 
depth of its gravitational potential. 


5.4. Intrinsic dynamics 

Our results demonstrate an anti-correlation between 
the residuals ALx^^^ and some intrinsic kinematic fac¬ 
tors such as rotation A, ellipticity e, and th e pararn eter 
for population A galaxies (see Figures A^ A8 and 
A9). Such an empirical pattern has been noted by previ¬ 
ous studies (e.g. Sarzi et al. 2010, 2013). Rotation could 
have taken away part of the kinetic energy of stellar mo- 
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Fig. 15.— Star formation rate as a function of cold gas mass, top-left: color code is the same as in Figure]^ top-right: “true lenticular” 
galaxies (brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black), bottom-right: 
field galaxies (blue) and galaxies in groups and clusters (magenta). 


tions and slowed down thermalization processes. Sarzi 
et al. (2010) found that the hot gas temperature of slow 
rotators met the stellar-thermalization kinetic energy ex¬ 
pectation while fast rotators fell short. Fast rotating sys¬ 
tems are able to expel hot gas to larger radii which leads 
to a more diffusely distributed halo (Negri et al. 2014). 
Meanwhile, it has been noted since the time of the Ein¬ 
stein that elliptical galaxies have larger Lx ILo^t than 
lenticular galaxies (Eskridge et al. 2005). Independent of 
the galaxy kinematical support, and apart from any evo¬ 
lutionary phase, flattening alone would reduce the bind¬ 
ing energy of galaxies of a given total mass. Yet two cor¬ 
relations, ALxgas" rotation and ALx^^^- ellipticity, are 
likely to be byproducts of one another. This is because 
rotation and flatness are inherently related. Rotational 


systems eventually become flatter along the major axis 
while flatter systems are able to possess higher rational 
support (Binney & Tremaine 1987). We find a signifi¬ 
cant anti-correlation between ALxgas ellipticity for 
slow-rotating galaxies, which favors the role of flattening. 
Sarzi et al. (2013) also found that a few Lx-deficient slow 
rotating galaxies appear to be flat. However, through 
two-dimensional simulations, Negri et al. (2013) demon¬ 
strated that flat galaxies and round galaxies can contain 
the same amount of hot gas if they have the same ro¬ 
tation, and a fast-rotating flat galaxy could also contain 
less hot gas than a slowly-rotating flat galaxy, which fa¬ 
vors the role of rotation. Interestingly, we find that both 
rotation and flattening are only relevant for population A 
galaxies. This mass-dependence has been noted by some 
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Fig. 16. — Crossing times as a function of times since gas removal 
for faint galaxies to reach average Lx /Lx through replenishment 
to new stellar-mass loss. The dashed black line indicates the equal¬ 
ity between these two time scales, [see the electronic edition of the 
journal for a color version of this figure.] 

simulation work (Negri et al. 2014, Posacki et al. 2013). 
More massive galaxies are more sensitive to flattening 
and rotation. 


5.5. Ram pressure stripping and stellar-mass loss 

In this work we investigated whether ram pres¬ 
sure stripping can account for the observed scatter in 
Lx^^/Lk- We compare the hot gas content per stellar 
mass (Lx^^^/Lk as well as ALx^^J and the strength of 
the instantaneous ram pressure relative to the gravita¬ 
tional restoring pressure (Pram/^grav), for a subsample 
of the 15 ETGs in the Virgo Cluster in Figure We 
observed an apparent correlation with a correlation co¬ 
efficient of p = —0.343 ± 0.073, reflecting the role of ram 
pressure stripping. 

Another important stripping process is due to the 
Kelvin-Helmholtz instability, which allows a shear force 
to be generated as intracluster gas passes by the outer 
edge of the galaxy; this shear force allows material to be 
continually stripped from the galaxy (Close et al. 2013). 
It is referred to as turbulence viscous stripping by some 
studies, as one of the transport processes (Nulsen 1982) 
which may cause substantial stripping of cluster galaxies 
at a stripping rate sometimes even greater than that due 
to instantaneous stripping alone. 

Galaxies move through the ICM in roughly elliptical 
orbits with one focus on the cluster center. As a galaxy 
approaches its pericenter, entering the inner (and denser) 
region of the cluster, it gains velocity towards the bot¬ 
tom of the gravitational potential well of the cluster. 
During this time it experiences stronger ram pressure 
stripping such that it has its hot gas stripped instan¬ 
taneously. When approaching its apocenter, reaching 
the larger radii of the cluster with a smaller velocity, 
it experiences weaker ram pressure stripping and stops 
losing gas. During this time, gas generated from new 


stellar-mass loss may be able to rebuild the atmosphere. 
While orbiting within the cluster, galaxies go through 
gas mass loss and gas mass replenishment phases alter¬ 
nately. NGC 4649 is an X-ray luminous galaxy in the 
Virgo Cluster that is not the center of a subgroup. We 
take NGC 4649 as a typical X-ray luminous early-type 
galaxy that has a benchmark Lx^^^/Lk- Twelve Virgo 
galaxies in this sample that have smaller /Lk than 

NGC 4649 may have lost some of their hot gas due to 
ram pressure. We estimate the time needed for these 
X-ray sub-luminous galaxies to gain enough gas to reach 
the average Mx^^^/Lk through stellar-mass loss. We cal¬ 
culate their crossing time using the radial vel ocity and 
the virial radius (Cr = ?'vir/'!^rad)- Figure 16 compares 


the time required for these X-ray sub-luminous galaxies 
to recover their gas mass with their crossing time. Ten 
out of 12 galaxies have a replenishment time shorter than 
their crossing time. These two time scales are compara¬ 
ble for the other two galaxies. This implies that most 
X-ray sub-luminous galaxies may regain their gas mass 
and become X-ray luminous galaxies at some point. 

6. SUMMARY 


This work presents our investigation of 42 FTGs in 
the ATLAS^d survey that have sufficiently deep Chan¬ 
dra observations to constrain their hot gas properties. 
In order to find the origin of the large scatter in the 
Lx — Lk relation of FTGs, we measure their hot gas 
content and relate it to other observables: stellar age, 
total galaxy mass, effective radius, depth of gravitational 
potential well, hot gas temperature, stellar velocity dis¬ 
persion, atomic gas mass, molecular gas mass, cold gas 
mass, star formation rate, angular momentum A, elliptic- 
ity e, A/i/e, environmental density, and relative strength 
of ram pressure (only for 15 Virgo galaxies). These sum¬ 
marize our main conclusions: 

• We obtain a scatter of up to a factor of 1000 in the 
Lx ^^~Lk relation for galaxies in our sample, consistent 
witE the results of previous studies. The scatter in the 
Mx^^^-Lk relation is slightly smaller than in the 

Lk relation. 

• We find two populations of early-type galaxies: 
galaxies with larger Lk are generally massive, slow- 
rotating, hot, big, and they tend to reside in high galaxy 
density environments, while galaxies with a smaller Lk 
tend to be low mass, fast-rotating, cool, and small galax¬ 
ies, which tend to reside in the low galaxy density envi¬ 
ronments. 

• The hot gas residuals (ALxg^^J are most strongly cor¬ 
related with the (rather low) star formation rates and hot 
gas temperatures; this indicates that stellar feedback in 
early-type galaxies heats up the gaseous component into 
the X-ray emitting phase and these galaxies are massive 
enough to keep such hot gas bound. Alternatively, both 
low level star formation and enhanced X-ray luminosities 
may be consequences of recent wet minor mergers. 

• We find that star formation rate increases with cold 
gas mass for low-mass, fast-rotating, and field galaxies, 
while massive, slow-rotating, and cluster galaxies tend 
to have constant SFR regardless of their cold gas mass. 
Various mechanisms may have interfered with the star 
formation process in the latter population. 

• Farly-type galaxies which contain significant 
amounts of cold gas, especially molecular gas, tend to 
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have more hot gas content. However, in our sample 
this trend only applies to “true-lenticular” galaxies, fast- 
rotating galaxies, low mass galaxies, and/or galaxies in 
low density environments. This could result from hot gas 
in the ISM mixing with the molecular gas. 

• Slower-rotating galaxies and rounder galaxies tend to 
have larger hot gas content. This trend only applies to 
slow-rotating galaxies, high mass galaxies, and/or galax¬ 
ies in groups and clusters. Flatness may play a bigger 
role than rotation. 

• While cold gas in cluster galaxies may be formed 

through hot gas cooling or stellar mass loss, the depo¬ 
sition time of the hot gas is longer than the stellar age 
of field galaxies. Cold gas in field galaxies is unlikely to 
be supplied by the radiative cooling of hot gas. Adding 
cold gas mass to the scatter in the Mx ^^-Lk 

relation for galaxies in this sample becomes even larger. 
This further suggests that the cold gas may be supplied 
by external mergers/accretions. However, the total gas 
mass Mgas increases almost linearly with Lx, which is ex¬ 
actly what we expect if stellar-mass loss is the primary 
source of all gaseous components. 

• For the 15 galaxies in our sample that are members of 
the Virgo Cluster, we established a correlation between 
the dehcit of their hot gas content and the strength of 
ram pressure, measured by the ratio of the instantaneous 
ram pressure to the gravitational restoring pressure. We 
demonstrate that stellar-mass loss is able to replenish 
stripped gas within the cluster crossing time for most 
galaxies. 

We conclude that the hot gas content per stellar light 


for a galaxy is the result of various mechanisms and 
feedback mechanisms. The situation varies for different 
galaxies with different evolutionary histories and envi¬ 
ronments. 
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Fig. A1.— Lxga.s/-^/c ^ function of stellar age (top-left), atomic gas mass (top-right), molecular gas mass (bottom-left), and cold gas 
mass (bottom-right). Color code is the same as in Figure]^ 
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Fig. A2.— as a function of star formation rate (top-left), hot gas temperature (top-right), total mass (bottom-left), and 

environmental galaxy density (bottom-right). Color code is the same as in Figure]^ 
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Fig. A3.— LxgasI ^ function of rotation {left) and ellipticity {right). Color code is the same as in FigurelJ^ 
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Fig. A4.— ALxgas ^ function of stellar age. top-left: color code is the same as in Figure top-right: “true lenticular” galaxies 
(brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black), bottom-right: field 
galaxies (blue) and galaxies in groups and clusters (magenta). 
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Fig. A5.— as a function of atomic gas mass, top-left color code is the same as in Figure|^ top-right: “true lenticular” galaxies 

(brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black), bottom-right: field 
galaxies (blue) and galaxies in groups and clusters (magenta). 
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Fig. A6.— ALxgas ^ function of cold gas mass, top-left: color code is the same as in Figure top-right: “true lenticular” galaxies 
(brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black), bottom-right: field 
galaxies (blue) and galaxies in groups and clusters (magenta). 
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galaxies (blue) and galaxies in groups and clusters (magenta). 
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Fig. A8.— ALxgas ^ function of angular momentum A. top-left color code is the same as in Figure [4| top-right: “true lenticular” 
galaxies (brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black), bottom-right: 
field galaxies (blue) and galaxies in groups and clusters (magenta). 
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Fig. A9.— ALxgas a function of ellipticity e. top-left: color code is the same as in Figure top-right: “true lenticular” galaxies 
(brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black), bottom-right: field 
galaxies (blue) and galaxies in groups and clusters (magenta). 
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Fig. AlO.— ALxgas ^ function of the environmental galaxy density, top-left: color code is the same as in Figure|^ top-right: “true 
lenticular” galaxies (brown) and “true elliptical” galaxies (green), bottom-left: high mass galaxies (red) and low mass galaxies (black). 
bottom-right: field galaxies (blue) and galaxies in groups and clusters (magenta). 
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B. STATISTICAL UNCERTAINTIES AND SYSTEMATIC TESTS 

The statistical uncertainties of these correlation coefficients were estimated through the Monte Carlo realizations of 
our measurements within their associated spectral fitting uncertainties (or uncertainties of data taken from literature, 
e.g. SFR). We listed such impact on each coefficient for the whole sample in Tables 4 and 5. The only significant 
impact is on relations with Tx due to its relatively large uncertainties of spectral fitting. The Tx-Lx Lk relation 
has p = 0.515 ± 0.184 and the Tjf-ALxgaa relation has p = 0.476 ± 0.186. 

The determination of the best-fit slope of the Lx^^^-Lk relation is very crucial to our study, since ALxg^a is very 
sensitive to this slope. In §4, we fit the Lx^^^-Lk relation to a power law log(Lxgaa) = 7l log(Lx) + B and obtain 
A = 2.3 ±0.3 and B = 14.6 ±3.3. The black solid line in Figure [^indicates the best-fit log(LxgaB) = 2.3 log(Lx) ± 14.6. 
We calculated ALxgaa again using A = 2.0 and A = 2.6 as indicated respectively by the dashed red line [log(Lxgaa) = 
2.0 log(L_R') ± 17-9] and dashed blue line [log(Lxgaa) = 2.6 log(Lx) ± 11-3] in Figure All Compared to the statistical 


uncertainties, the change of A by ±0.3 has little impact on our results as listed in Tables 4 and 5. Boroson et al. (2011) 
determined a best-fit of A = 2.6 ± 0.4 using a smaller sample of galaxies. We also performed the fit of the Lx^^^-Lk 
relation ignoring those poorly constrained datasets (8 extremely faint galaxies with their Axg^a consistent with zero) 
and we obtained A = 2.1 ± 0.3. Both these cases are within our variations. 

We also performed the maximum likelihood procedure by Kelly (2007) to quantify the Lx^^^-Lk relation. It is a 
Bayesian method that includes both individual measurement error and population intrinsic scatter within the linear 
regression. The best-fit determined by its linmix_err regression is log(Lxg,,J = (2.2 ± 0.3) log(Lx) ± (15.8 ± 3.3), in 
fully consistent with the best-fit we obtained. It also falls into the range of variations we examined. We also listed the 
slope of the Lx^^^-Lk relation determined by liiimix_err for each subgroup in Table 6, in agreement with the results 
we obtained. 

The Lx-Lopt relation obtained by O’Sullivan et al. (2001) follows a broken power law (bknpower) (see Figure!^. 
That may be a result of them including the LMXB components and biasing high the luminosities of faint galaxies due 
to the limited resolutions of ROSAT. We also fit the Lx^^^-Lk relation of galaxies in our sample obtained with high 
quality Chandra observations to a bknpower and obtained a best-fit of 


lotrlLv 1 - I + 

log(±Ag.J - I 3^iiog(Lx) ± 


27.0 (log(Lx) < 10.7) 
5.4 (log(Lx) > 10.7) 


The best-fits of the two slopes are 1.1±0.1 and 3.1±0.5. They are only different by 1.8 (t. Furthermore, we compared 
the standard deviations of the distributions of ALxgas obtained by using bknpower and single power law respectively: 
0.690 ± 0.07 and 0.737 ± 0.109. They are not significantly different from each other. We also performed ftest and we 
do not find bknpower improve the fit significantly. Therefore, we do not think it is statistically necessary to employ 
bknpower. Still, we tested our results using a bknpower as the best-fit of the Lx^^^-Lk relation as listed in Table 5. 
Its impact does not change our results qualitatively. The “V-shaped” relation between ALx^.^^ and Mtot (’’e, Aftot/ce, 
and cr) diminished (Figure A12 left), while ALx^^^ is still most strongly correlated with star formation rate and hot 
gas temperature with even larger correlation coefficients (see Figure [Al2| -right) . 

In addition to the Spearman correlation, we also use Kendall’s tau correlatiorj^ to quantify the relations between 
various factors. None of these results are qualitatively different from using the Spearman correlation. The tau 
correlation coefficient of each relation is listed in Tables 4 and 5. 


\protect ittp://cran.r-project.org/web/packages/NADA/NADA.pdf 
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Fig. All.— Gaseous X-ray luminosities as a function of iC-band luminosities for all galaxies in this work. The black solid line is the 
best-fit relation ( log(Lxga.s) = 2.31og(Lx) -h 14.6). Dashes black line is the best-fit Lxgag-T/^ relation to broken power law. 

Dashed red line indicates log(Txgas) ~ 2.01og(Lx) + i-7.9 and dashed blue line indicates log(Txgas) “ 2.61og(Lx) + 11.3. Color code is 
the same as in Figure]^ 




Fig. A12.— ALx, 
best-fit of the Lv 


as a function of total mass (left) and star formation rate (right) when using a broken power law function as the 
-Lx relation. Color code is the same as in Figure 
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